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ABSTRACT 

Currently,  topographic  products  such  as  contour,  slcpe,  and 
aspect  are  commonly  derived  from,  computer  processing  of 
digital  grid-point  elevation  models  (DEM).  OEM's  can  be 
laboriously  created  by  a  user  from  an  analytical  plotter 
equipped  with  profiling  capability.  Alternatively,  a  user 
can  acquire  e  DEM  produced  by  an  external  source,  bi.t  then 
be  constrained  by  the  production  resolution  or  accuracy  which 
night  not  be  suitable  for  the  user's  application. 

Recently,  the  integration  of  stereoplottcrs  with  geographic 
information  systems  has  resulted  in  limited  terrt in  data 
production  in  X,Y,Z  coordinates  (3-d inens ions  1  cigital 
cartegraphic  data).  In  those  eases  where  detailed  3-dim >nsional 
date  exists,  building  a  DEM  from  this  database  can  stve  labor 
and  preserve  user-defined  resolution  and  accuracy  criteria. 

In  this  study,  previously  digitised  3-dinen s ions  1  terrain 
data  are  combined  and  converted  to  a  3  arc  second  (user ‘defined) 
cell  matrix.  Multiple  elevation  data  points  falling  into 
the  same  call  were  averaged  to  provide  a  single  e  Ovation 
value  aad  these  cells  without  an  elevational  value  were 
flagged.  The  existiag  elevation  date  could  then  be  lisplayed 
via  stereo  superposition  over  a  photographic  model  established 
ea  aa  analytical  platter  sad  those  flagged  cells  <ould  be 
profiled.;  This  paper  presents  the  method  involved  in  building 
a  DEN  from  3-dimeasioaal  cartographic  data  and  weigLs  the 
benefits  and  shortfalls  of  this  method. 


INTRODUCTION 


To  date,  spatial  data  capture  and  representation  of  a 
3-dimensional  physical  world  is  accomplished  primarily  through 
a  2-d imens iona 1  (X  and  Y)  data  capture  process.  Elevation, 
(Z),  has  largely  been  treated  as  a  separate  entity  unrelated 
to  physical  features  on  the  topography.  Recent  attempts  to 
compensate  for  this  deficiency  have  focused  on  the  merging 
of  planimetric  data  with  digital  grid  point  elevation  data. 
The  merging  of  Landsat  data  with  Defense  Mapping  Agency  (DMA) 
digital  terrain  elevation  data  ( DTED )  (Tanaka,  1979)  or  DMA 
digital  feature  analysis  data  (DFAD)  with  DTED  (Federhen, 
1984)  illustrates  recognition  of  this  problem. 

OEMs  provide  a  relatively  compact  and  efficient  machine-processable 
format  for  representing  physical  topography.  Contours,  slope, 
and  aspect  are  commonly  derived  from  computer  processing  of 
DIMe.  A  user  requiring  digital  terrain  elevation  data  must 
either  acquire  a  DEM  produced  by  an  external  source  or  create 
this  product.  Frequently,  a  DEM  is  not  available  from  an 
external  source  for  a  given  area.  Futhernore,  if  one  does 
exist  it  may  not  be  at  the  required  resolution  or  accuracy 
needed  for  a  specific  application.  Alternatively,  a  DEM  may 
be  laboriously  created  from  an  analytical  plotter  equipped 
with  a  profiling  capability.  Errors  are  possible  and  can  be 
introduced  in  this  process  due  to  the  redundancy  of  the 
operator’s  task  and  to  hardware  plotter  problems  encountered 
during  profiling. 

■oedat  merging  of  stereoplotters  with  geographic  information 
eyatoma  has  resulted  in  limited  terrain  data  productions  in 
X,  Y,  Z  coordinates.  This  permits  direct  data  entry  by  an 
analyst  and  storage  of  this  data  in  a  coaputer-processable 
format.  Growth  of  this  technology  and  production  of  3-D 
digital  data  will  intensify  due  to  market  demand  from  such 
diverse  arose  as  robotic  vehicle  navigation,  computer  image 
generation,  and  terrain  modeling.  In  those  cases  where 
3-dimemsioasl  data  exists,  building  a  DEM  from  this  database 
can  save  labor  and  preserve  user-defined  resolution  or  accuracy. 

In  this  study,  a  detailed  3-D  terrain  analysis  database  was 
used  as  the  basis  for  the  production  of  a  digital  grid  point 
elevation  database.  This  paper  briefly  describes  the  hardware, 
software,  and  methods  used  to  generate  this  DEM. 


HARDWARE  and  SOFTWARE 


Efforts  at  the  USAETL,  through  a  program  in  computer-assisted 
photo  interpretation  reseach  (CAP1R),  have  matched  an  analytical 
platter  equipped  with  stereo  superposition  graphic*  with  a 
geographie  information  system  (CIS)  to  provide  the  mechanism 
for  3-dimensional  data  capture  and  management.  The  analytical 
plotter  used  daring  this  study  is  APFS-IP.  It  is  a  medium 
aesnrasy  plotter  rated  at  tea  micrometers  sfter  correction 
nsiag  at  affine  transf ornat ion  (Crsves,  1982).  Stereo 
superposition  graphic  CRT's  are  mounted  in  the  rear  housing 
of  the  plotter.  The  graphic  CRTs  display  digitised  spatial 
data  which  is  introduced  through  bean  splitters  direetly  into 


t he  optical  path  of  the  plotter  stereoscope.  This  is 'viewed 
by  the  analyst  3-diaensionally  superimposed  over  the 
photograaaetric  aodel.  These  graphics  are  updated  ten  times 
per  second  so  that  the  displayed  data  tracks  the  imagery  as 
the  analyst  traverses  the  stereo  aodel.  The  APPS-IV  has  beer, 
enhanced  with  a  profiling  capability  which  allows  collection 
of  digital  elevation  data  along  evenly  spaced  profiles. 
Editing  a  portion  of  the  digital  elevition  matrix  can  be 
perforaed  within  the  profiling  system  without  having  to 
reconstruct  the  entire  DEM. 

The  geographic  information  system  used  for  data  capture  during 
this  effort  was  the  Analytical  Mapping  System  (AMS).  It 
allows  one  to  aaro-tr iangulate ,  digitize,  edit,  verify,  and 
then  database  the  results.  The  digitization  is  in  arc/node 
format  with  attribute  entry  and  the  data  is  stored  in  geographic 
coordinates . 


3-D  TERRAIN  DATABASE 

In  this  study,  3-diaens ions  1  digital  cartographic  data  was 
taken  froa  a  terrain  analysis  of  the  Fort  Belvoir  area 
conducted  in  order  to  test  and  demonstrate  CAPIR  capabilities 
developed  at  DSA8TL  (Edwards, 1983) .  NASA  0-2  color  infrared 
photography  (scale  1:130,000)  was  the  priaary  data  source 
far  the  terrain  analysis.  Two  stereo  pairs  were  triangulated 
using  field-checked  and  photo-identified  first  and  second 
order  survey  points.  Initially,  landfora,  surface  drainage, 
sails,  and  landcover  were  individually  interpreted  and 
digitised  ae  shown  in  Figure  1.  A  transportation  classification 
wee  later  started  in  conjuction  with  this  study,  but  currently 
ramaiae  partially  completed.  The  actual  terrain  category 
boundaries  in  each  classification  were  :  nterpreted  during 
ea-liae  atereo  digitisation.  The  graphi:  feedback  provided 
by  stereo  aaperposit ion  of  the  cartograpnic  data  permitted 
visual  verification  of  the  geographic  accuracy  of  the  data 
during  digitisation. 

The  terrain  aaalyais  was  limited  to  an  area  covered  by  the 
8*8*  Geological  Survey  Fort  Belvoir  7  1/2  ainute  quadrangle. 
A  3  are  eecoad  OEM  for  this  area  contains  22,300  elevation 
points.  Thin  could  be  laboriously  created  using  the  profiling 
capability  of  the  analytical  plotter.  I-i  this  investigation, 
the  existing  elevation  data  within  the  five  listed  terrain 
analyaia  files  were  combined  to  form  a  partial  DEM.  Profiling 
could  ebon  used  to  complete  the  DEM  matrix. 


DIM  SYNTHESIS 

The  following  eoction  details  the  steps  in* olved  in  reformatting 
the  cartographic  data  fron  a  vector  foraat  (as  produced  by 
AMS)  into  a  raster  format,  with  the  grii  spacing  of  the  DEM 
defined  by  the  user.  This  conversion  was  accomplished  by  a 
dombiaation  of  Fortran  programs  and  lort/Merge  routines 


The  data  onperted  fron  AMS  is  changed  fton  an  arc-node  format 
to  a  polygon  format  complete  with  a  polygon  attribute  header 

_ 


Drainage 


tag.  In  order  to  use  this  data,  the  header  information  had 
to  be  stripped  off  the  polygon  data  files  leaving  only  poinc 
date.  The  number  of  3-dimensional  data  points  in  each  of 
the  five  terrain  analysis  files  varied  widely  and  are  listed 
below. 


Terra:.n  Analysis  File 
Landf orm 
Drainage 
Soil 

Landeover 
Transportat ion 


Humber  of  Points 
37  15 
15325 
7543 
15665 
776  1 


In  order  to  simplify  the  programming,  all  of  the  points  were 
merged  into  one  file  of  50009  points.  The  data  was  then 
sorted  by  latitude  and  the  minimum  and  maximum  values  stored 
to  provide  the  northern  and  southern  boundaries  of  the  DEM. 
■eat,  the  %-d inens iona 1  values  were  sorted  by  longitude, 
placed  into  3  arc  second  divisions,  and  then  resorted  by 
latitude  within  each  of  these  longitudinal  divisions.  The 
longitude  vi.lues  were  sorted  in  descending  order  because 
longitudinal  values  decrease  from  west  to  east  in  the  western 
hemisphere.  This  sorting  provided  an  ordered  grid  of  X,Y,Z 
data  sorted  oy  latitude  and  longitude. 

After  sorting,  multiple  entries  of  the  same  X,Y,Z  data  points 
were  evident  This  is  a  result  of  the  AMS  software  transforming 
arc-node  data  into  polygon  data.  Points  which  are  members 
of  boundary  lines  between  2  or  more  polygons  are  entered  for 
each  of  the  individual  polygons.  For  example, 


Figure  2. 

in  Figure  2  the  latitude,  longitude  and  elevation  for  points 
1,2,3,  and  A  will  appear  for  both  polygons  L  and  M.  The  data 
for  point  3  appears  in  polygons  L,  M  and  R.  The  redundant 
data  values  were  deleted  for  two  reasons.  First,  thi*  reduced 
the  overall  amount  of  data  to  be  handled  and  thus  decreased 
the  computer  processing  tine.  Second,  in  a  later  program 
the  elevations  which  are  included  in  each  DEM  grid  cell  (3 
are  second)  are  averaged  to  compute  the  elevation  for  that 
cell.  Multiple  entries  of  the  same  value  will  bias  this 
resale,  j ta  the  Fort  Selvoir  test  case,  ehe  number  of  data 
Mints  Ms  reuueed  from  50009  poines  to  24157  points.  This 
reduced  dats  was  then  combined  to  create  a  partial  DEM. 
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Within  this  partial  DEM,  if  a  grid  cell  had  2  or  more  elevation 
values,  then  the  cell  was  assigned  the  average  of  the.* 
values.  If  no  data  was  available  then  that  particular  cell 
was  tagged  with  the  value  99999.  This  formed  a  matrix  of 
values,  a  partial  DEM. 

The  areas  within  this  DEM  containing  no  elevation  value  :an 
then  be  filled  in  using  the  profiling  capability  of  the 
APPS-IV.  The  partial  DEM  could  be  displayed  using  the  graphic 
superposition  CRTs  mounted  in  the  rear  of  the  APPS-IV.  The 
cells  containing  no  elevational  value  (99999)  would  not  be 
displayed  over  the  photogranmetric  model  established  on  the 
plotter.  Then,  using  the  profiling  capability  of  the  analytical 
plotter,  the  areas  without  elevation  data  could  be  bounded 
and  profiled. 


DISCUSSION 

This  method  of  creating  a  DEM  from  3-dimensional  cartographic 
data  has  several  advantages.  First,  it  exploits  available, 
existing  data.  Second,  the  time  and  effort  required  to  create 
a  DIM  is  minimized  since  only  a  portion  of  the  matrix  must 
be  compiled  through  profiling.  Finally,  the  accuracy  of  the 
DBM  may  be  improved  since  the  averaging  of  the  elevation 
values  within  a  cell  reinforces  that  final  value. 

Several  limitations  are  inherent  in  this  method.  First,  this 
process  requires  3-dimensional  cartographic  data.  Currently, 
Z,Y,Z  cartographic  data  is  not  commonly  available.  However, 
as  mentioned  earlier,  production  of  this  data  is  expected  to 
become  more  widespread  as  additional  systems  linking  an 
analytical  plotter  with  a  CIS  become  available.  Second,  the 
benefits  of  using  this  method  increase  with  the  density  of 
the  cartographic  data  set.  Clearly,  the  more  detail,  the 
less  data  needs  to  be  gathered  through  profiling.  Sparse 
data  sets  accrue  negligible  benefits  from  this  method, 
finally,  the  averaging  of  elevation  values  within  a  cell  has 
a  smoothing  effect  on  the  DEM.  Thus,  extreme  values  are 
filtered  out  in  these  cases. 

A  fairly  simple  method  of  DEM  synthesis  has  been  preset  ted 
in  this  paper.  In  the  future,  more  sophisticated  techniques 
will  enhance  this  procedure.  Artificial  intelligence  rule-btsed 
elassifiers  which  would  relate  elevational  values  bases  on 
select  terrain  category  types  such  as  water,  plains,  or 
alluvial  fans  would  further  reduce  the  DEM  synthesis  effort. 
Such  rule-based  classifiers  would  be  handled  by  an  e;  pert 
system. 

* 

In  conclusion,  s  variety  of  methods  and  tools  will  be  reeded 
to  meet  the  need  for  terrain  infornation.  In  the  fu:ure, 
advances  will  he  made  in  cartographic  data  capture,  managenent, 
nnd  exploitation.  This  techniquo  of  using  existing  3-dimen.ional 
data  to  events  a  DIM,  has  potential  to  be  one  of  these  tools. 
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